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Abstract We investigate the possibility that prebiotic homochirality can be achieved 
exclusively through chiral-selective reaction rate parameters without any other 
explicit mechanism for chiral bias. Specifically, we examine an open network of 
polymerization reactions, where the reaction rates can have chiral-selective values. 
The reactions are neither autocatalytic nor do they contain explicit enantiomeric 
cross-inhibition terms. We are thus investigating how rare a set of chiral-selective 
reaction rates needs to be in order to generate a reasonable amount of chiral bias. 
We quantify our results adopting a statistical approach: varying both the mean 
value and the rms dispersion of the relevant reaction rates, we show that moderate 
to high levels of chiral excess can be achieved with fairly small chiral bias, below 
10%. Considering the various unknowns related to prebiotic chemical networks in 
early Earth and the dependence of reaction rates to environmental properties such 
as temperature and pressure variations, we argue that homochirality could have 
been achieved from moderate amounts of chiral selectivity in the reaction rates. 

Keywords Prebiotic Chemistry • Chirality • Exoplanets 



1 Introduction 



The origin of life on Earth is a central scientific question of our time | Sullivan 
[and B aross} |2007| |Orgel| |2010|). A possible way of organizing its many facets is 



to separate it into three basic aspects, all still unknown: when did it happen; 
where did it happen; and how did it happen. The "when" question is, to a certain 
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extent, the best understood. We can put a lower limit on the origin of life on Earth 
at about 3.5 billion years ago (BYA) (Mojzsis et al. 19961. Possible evidence of 
earlier life, dated at 3.8 BYA, remains controversial (van Zuilen, Lepland, and 



Arrhenius 2002 I. Of course, it is also possible that life emerged more than once 



before this time, although if it did it was probably made extinct quickly due to the 
unstable conditions of a crustless primeval Earth. Another intriguing possibility 



is that of a shadow biosphere, which has hitherto remained undetected (Benner, 
Ricardo, and Carrigan} |2004| |Davies et al.[ |2004[ ). At about 3.9 BYA, the Late 



Heavy Bombardment caused by the displacement of outer planets, constitutes a 



serious obstacle against any long-term survivability of earlier life forms ( Corners 



et al. 



2005). 



The "where" question has many possible answers, ranging from shallow wa- 



ter pools and hot water vents to clay and mineral surfaces ( Corliss, Baross, and 



Hoffman 1981 ). It is also speculated that life may have originated in outer space. 



possibly Mars (|McKay 20101, and was transplanted to Earth via meteoritic im- 
Crick and Orgel! 



pacts (Davies 



1996 



1973 



Sullivan and Baross 2007 1 . If not life 



itself, at least some of its basic organic ingredients may indeed have come from 



outer space, possibly seeding early life here ( Chyba and Sagan 

2008J . 



Pizzarello 



1997 Martins et al. 



1992 



Cronin and 



The "how" question remains quite open and it may never be answered precisely, 
as it depends on the details of a prebiotic biochemistry we may never completely 
know. However, the uniformity in the biochemistry of all living organisms strongly 
suggests that all life descended from a last universal common ancestor (LUC A), 
that lived 3.2-3.8 BYA and had a level of complexity similar to that of a simple 
modern bacterium ( Orgel[ 2010). To "explain" the origin of life requires us to ex- 
plain how abiotic organic compounds combined to generate the LUCA, at present 
a major challenge. 



A key aspect of all terrestrial life is its biomolecular homochirality (Bonner 
1995 Fitz et ah] 2007[ ): the vast majority of amino acids comprising living creatures 



are levorotatory (L), while the majority of sugars are dextrorotatory (D). Civen 
the link of present and past life to the LUCA, we should thus expect that the 
LUCA displayed the same chiral signature as modern life. To explain the origin 
of life we must also therefore explain the origin of homochirality. Although it is 
possible that biomolecular homochirality evolved during, or even after, prebiotic 
chemistry reached a more complex stage as, say, when the first nucleic acids were 
being formed ( jCilbert 1986') , it is more likely that homochirality is a precondition 
to life. Accordingly, we will assume here that chiral selectivity played a key role as 
prebiotic biochemistry began to evolve, setting the conditions for the first living 
forms to emerge. Within this framework, to understand the "how" question we 
must first understand the homochiral question: essentially, we must understand 
how the first homochiral polymers formed under prebiotic conditions. 

Over the past decades, several explanations have been offered to account for 
the homochirality of living organisms (Bonner) ,1995 Fitz et al. 2007 1 . All of them 
assume that the essential reactions leading to homochirality are highly nonlinear 
and unstable to growth, so that a small excess of one handedness can be efficiently 
amplified to full homochirality. There are thus two basic aspects of the homochiral 
question: first, the cause for the initial excess of one handedness over the other; 
and second, the amplification mechanism. 



Chiral Polymerization in Open Systems From Chiral-Selective Reaction Rates 



3 



Causes for the initial chiral excess can be classified as endogenous or exoge- 
nous. By endogenous causes we mean those that are intrinsic to the physics and/or 
chemistry of the reactor pool. For example, parity symmetry breaking in the weak 
nuclear interactions has been proposed as a possible source of the initial homochi- 



ral excess ( Yamagata[ 



1966 



Salam 



1991 



Kondepudi and Nelson 1985 1 , although 



the numbers obtained in quantum computations using perturbed Hamiltonians are 
exceedingly small (Bakasov, Ha, and Quack 1998). In contrast, exogenous sources 
are caused by outside influences. Examples of exogenous sources of initial chiral 



bias include circularly-polarized UV light from active star formation regions ( Bai- 



ley 2001 1, and deposits of chiral compounds by meteoritic bombardment (Cronin 



and Pizzarello 



1997 



Glavin and Dworkin 2009). Another possibility of an ex- 



ogenous influence is that chirality is induced by a clay template, such as in recent 



20111, which exhibit 



results for n-Propyl NH3CI vermiculite clay gels (Fraser et al. 
chiroselective absorption of amino acids. 

At this point, the debate is still open as to whether any of these possible sources 
of chiral bias was present at prebiotic times and, if it was, whether it was efficient 
or necessary. Gleiser, Walker, and Thorarinson have recently argued that even 
if an initial chiral excess was produced by any of these (or other) mechanisms, 
environmental disturbances could have erased it, restoring either achiral condi- 
tions or switching chirality in the opposite direction: chiral selection might be the 
result of a series of stochastic reversals prompted by large temperature and/or 



density fluctuations, or other possible environmental effects (Gleiser, Thorarinson, 
and Walker 20081. Similar conclusions were reached recently by Hochberg ( 2009[ ). 



Accordingly, the question of biological homochirality cannot be separated from 
the details of planetary history. The same holds for any other planetary platform 
displaying stereoselective chemistry. 

The second basic aspect of the homochiral question focuses on the details 
of the amplification mechanism. Ever since the seminal work by Frank (Frank 



1953), most models attempt to describe how a small chiral bias can be efficiently 



amplified toward homochirality. Such models rely on a combination of two proper- 
ties: the reactions must be autocatalytic, thus guaranteeing the nonlinear behav- 
ior needed for amplification, and they must display enantiomeric cross-inhibition, 
that is, enantiomers of opposite chirality can combine, suppressing the growth of 



homochiral chains (Sandars 



2003 



Wattis and Coveney 



2005). As can be seen in the Sandars model (Sandars 20031, this inhibition may 



2005 



Saito and Hyuga 



occur by suppressing the enzymatic activity of homochiral chains in the reactor 
pool. Within this framework, one must still appeal to a mechanism for generating 
an initial excess. Traditional approaches along these lines have held that the initial 
symmetry breaking is endogenous, occurring at the level of small stochastic fiuc- 
tuations in the concentrations of left— and right— handed monomers, which must 
then be amplified by the combined effects of autocatalysis and cross-inhibition 



2003 



Gleiser and Walker 



to yield an appreciable enantiomeric excess (Sandars 
[2008'). However, fiuctuations between the concentrations of left- and right- handed 
enantiomers may be exceedingly small in a large reactor pool, thereby requiring 
strong autocatalytic feedback for amplification to occur, if at all. Given that au- 
tocatalytic feedback leading to significant chiral selection has only been observed 
in a few abiotic systems ( Saito and Hyuga} 2005), and has not yet been observed 
in a prebiotically relevant scenario, it is prudent to consider other possibilities for 
the onset of chiral selection in a prebiotic environment. 
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In this work we therefore explore an alternative scenario for the generation 
of a chiral excess in open polymerization systems without invoking autocatalysis 
or enantiomeric cross-inhibition. Instead, we propose that reasonable amounts of 
chiral excess can be achieved solely through chiral-selective reaction rates. In or- 
der to sample a large number of possible values for the reaction rates, we present 
a statistical analysis whereby each run represents a potential prebiotic chemistry 
with fixed values of the reaction rates. For each prebiotic scenario the rates are 
picked from a sample of Gaussian-distributed values. Although a large statisti- 
cal sampling of the systems under study would on average be racemic, we show 
that moderately small chiral-selective variations in reaction rates may lead to 
significant chiral excess in a steady-state. Within a prebiotic scenario on early 
Earth, we might imagine that such chiral-selective fluctuations in the reaction 
rates were prompted by environmental disturbances or by still unknown details of 
the prebiotic chemistry. We propose a few possibilities, although the strength of 
our statistical analysis is its independence from the specific mechanism generating 
the chiral excess in the reaction rates. Thus, we will explore the possibility that 
terrestrial biomolecular homochirality might have arisen without the need for any 
of the conventional mechanisms for chiral biasing or amplification. For illustrative 
values of the relevant reaction rates, we will place lower bounds on the amplitude 
of fiuctuating reaction rates needed for substantial chiral bias. 



2 Chiral Polymerization Model 



Consider a polymerization reaction network where activated levorotatory and dex- 
trorotatory monomers {Ll and Dl, respectively), can chain up to generate longer 
homochiral or heterochiral molecules. (Here, "activated" refers to these monomers 
being highly reactive due to energy input into the reactor pool, as in the APED 
model of Plasson et al, proposed as an alternative to Frank- like models (Plasson 



et al. 2004 1). The model reactions include deactivation of activated monomers. 



and the polymerization reactions, 

Ll + Li Li+i 
Ll + Dj °^ Mij 
Ll+M,, "A M,+ij 



Dl H Dl 

Dl + D, ^ 
Dl + L,"-^ M.i 
Dl + M„- ^ M,,+i 



(1) 

(2) 
(3) 
(4) 



where /il(_d) is the deactivation rate for L(Z?)-monomers and a^j^j) is the poly- 
merization rate for adding activated L(D) -monomers to a growing chain. Here Li 
and Dj denote homochiral polymers of length i and j respectively, and Mij de- 
notes polymers of mixed chirality consisting of i L-monomers and j D-monomers. 
Note that this reaction network is substantially different from the usual Frank-like 
model for chiral symmetry breaking ( Frank[ 1953 Sandars 20031: there is no au- 
tocatalysis (either explicit or through enzymatic activity of homochiral chains) or 
enantiomeric cross-inhibition present. We note that an alternative copolymeriza- 
tion model with two chemically distinct monomers has been recently proposed by 



Wattis and Coveney iWattis and Coveney 20071 
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What kinds of situations could lead to chiral-selective values for the reaction 
rates? It is well-known that the temperature dependence of reaction rates is de- 
termined by the Arrhenius factor, k ^ exp[—Ea/RT], where Ea is the activation 
energy related to departures from equilibrium populations of molecules. Consider, 
for instance, a situation where levorotatory and dextrorotatory rates have differ- 
ent temperature dependencies, parameterized by chiral-specific activation barriers 
as ki^(^£)-^ = exp[— _E^(£))(T)/_RT]. For example, we could model this dependence 
as -E_L(D) = — c^(D)T^ + dL{D}j where the four constants C£(d) s-nd dj^^jj-^ are 
positive. If, say, cl > cd and d^ = do, we see that for all temperatures (for 
T < = dn/cD) dextrorotatory reactivity is enhanced. A different realization 
would have correspondingly different constants or a different dependence of the 
activation energy with temperature. Each choice would correspond to an inde- 
pendent prebiotic scenario. (We should consider this model only as a suggestive 
mechanism whereby chiral selectivity could occur through the reaction rates.) 

Given that the details of prebiotic chemistry remain elusive, here we are con- 
sidering that, in principle, reaction rates may be chiral-selective. However, instead 
of assigning ad hoc values for the rates of opposite chirality, we will conduct a 
statistical analysis whereby, for different numerical experiments-that is, different 
realizations of prebiotic scenarios-the values for L and D reactions rates will be 
allowed to randomly fluctuate about a given mean. In other words, each time we 
solve the coupled nonlinear ordinary differential equations describing the reactions, 
we will pick a set of random, Gaussian-distributed values for the four reactions 
rates /il(_d) and 0^(15). This heuristic approach will allow us to investigate the 
range in the values of the reaction rates that generates a substantial departure 
from racemic results. As we will see, modest values for the fluctuation range for 
L and D reaction rates will suffice to generate substantial chiral bias. Although 
we do not model any specific chemical system, we expect that our general re- 
sults will be applicable to a range of experiments probing the rate-dependence of 
polymerization models and its impact on chiral symmetry breaking. 

In practice, we will consider pairs of randomly-distributed reaction rates (/iz,, /id) 
and (ttL, an) taken from a Gaussian distribution with a racemic mean 



where ho and ao are the rms values, and 5xx' is the Kronecker delta. The angled 

brackets mean an average over A'^ experiments (an ensemble average), with A'^ 
large. For example, the normalized probability P{hL) that the L-hydrolysis rate 
will obtain a value Kl is 



{Hl) = {ho) = h 
{ai) = {ao) = a, 



(5) 



and a two-point correlation 



{o-L{D)0-L'(D')) = (H)5lL'{DD'), 



(6) 



PihL) 



1 



exp[-(/ii - h)^ /2hl]. 



(7) 



For the above system, the rate equations for the various concentrations are 
written as (from here on, quantities in capitals denote the concentrations of the 
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different reactants): 



dpi 
dt 

dLl 

dt 
dPi 

dt 
dLj 

liF 

dt 
dMij 
dt 



S-hLLl-aLLl K]Li + ^D,+^^Mi,- -dLl (8) 

\i=l J = l i=lj = l / 

(CO OO OO CXJ \ 

^L, + ^Dj+^^My -dDl (9) 
i=i i=i i=ij=i ) 

hi^Ll — ulLiLI — aoLiDl — dLi (10) 

hoDl - ulDiLI - aoDiDl - dDi (11) 

GLLlLi-i - ttLLiLl - aoLiDl - dLi {i > 2) (12) 

aoDlDj-i - aLDjLl - aoDjDl - dDi {j > 2) (13) 
aL{LlM,-ij - MijLl) + aoiDlMij-i - MtjDl) - dM^ > 1014) 



Note that for the open systems studied here we also include the additional 
rates S describing the permeation of activated Ll and Dl monomers into the 
system, and d describing the loss of all molecular species from the system. For 
closed systems S = d = 0. 

We can simplify the reactions above by defining 

OO OO OO OO 

i=2 j=2 i=lj = l 

and summing the reaction equations over i and j. Then the reaction-network, in 
the limit of infinite-length polymers, reduces to: 



dLl 
dt 
dpi 
dt 
dLl 
liF 
dPi 
"dT 
dC 
dt 
dp 
dt 
dM 
dt 



S-HlLI - ulLI (Li+C + Di+V + M)- dLl; (16) 

S - hoDl - aoDl {Li + C + Di +V + M) - dDl; (17) 

HlLI — ulLiLI — aoLiDl — dLi; (18) 

hoDl - ulDiLI - aoDiDl - dDi; (19) 

ulLILi - {aoDl + d)C; (20) 

aoDlDi - (ulLI + d)V; (21) 

ulLI {Dl +V) + aoDl (Li + C) - dM. (22) 



The quantity of interest is the enantiomeric excess [ee(t)], that is, the net 
chirality, which we define as 

... ^ im) + Li{t) + m] - [DKt) + Di{t) + v{t)] 

^> [Ll{t) + Li{t)+C{t)\ + [Dl{t) + Di{t)+V{t)y ^ ' 



Chiral Polymerization in Open Systems From Chiral-Selective Reaction Rates 



7 



Note that our definition of the enantiomeric excess implicitly assumes that we 
are only interested in pure homochiral polymers. Given initial values for the con- 
centrations, the equations describing their evolution can be solved and the net 
chirality computed as a function of time. We will choose the simplest possible 
initial conditions, setting all initial concentrations to zero, with the exception of 
L^(0) and Dl{0), which we will set equal to each other so that the system starts 
with no chiral bias and no long polymer chains. We will also pick a constant value 
for the source for activated monomers, S{t) = S. With this choice, we will be able 
to investigate under which conditions, that is, under which values for the reaction 
rates, the system engenders a net chirality as it evolves to a steady-state wherein 
all concentrations achieve a dynamic equilibrium so that their time derivatives 
vanish. We note that results were not very sensitive to the choice of S. 

Before we show solutions to the system of equations, we eliminate the loss rate 
from the equations, dividing them on both sides by d. Since d has dimensions 
of inverse time, this is equivalent to introducing the dimensionless time variable 



T = td such that we can write the equations as, 

^ = a - I3lLI - olLI {Li+C + Di+V + M)- Ll; (24) 

UT 

^=a- PdDI - aoDl [Li + C + Di + V + M) - Dl; (25) 

UT 

= PlLI - olLiLI - aoLiDl - Li; (26) 

ar 

^ = I3dDI - aLDiLl - aoDiDl - Dl- (27) 
dr 

^ = aLLlLi - {aoDl + 1)C; (28) 
ar 

^ = aoDlDi - {aLLl + 1)V; (29) 
ar 

^ = aLLl{Di+V) + aDDl{Li+C)-M, (30) 

where we have introduced a = S/d, aL(^E)) = o,L{D)/d, and j3i^(^D) = hj^i^Q-j/d. 
With this parametrization uq = ao/d, fio = ho/d, and a = a/d and j3 = h/d 



characterize the random distributions of reaction rates outlined in eqs. [5] and ID 
Thus, all reaction rates and dynamical time-scales are expressed in terms of the 
disappearance rate d. In particular, a gives the ratio per unit volume of entrance 
versus exiting of reacting materials to and from the reactor pool. 

3 Results 

In Figure [l] we show the behavior of the time evolution of enantiomeric excess for 
several illustrative examples with a = 200 and different choices for the reaction 
rates ql, qd, Pl and Pd- On the left, we show results with = Pd = 100 fixed, 
and varying ratios of ul/oid- Shown are the results for ul = 5.0, an = 10.0 
{ciL/ctD = |) in blue, = 5.0, ao = 15.0 {aL/ao = i) in black, and Ul = 5.0, 
an = 20.0 (ai^jau = \) in red. The system quickly reaches a steady state with net 
chiral excess eCss — 0.24, ecss — 0.38, and ecss — 0.48, respectively (corresponding 
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to 24%, 38%, and 48% enantiomeric excess). On the right, we show results holding 
oiL = cxd = 10.0 fixed, and varying the ratio of Pl/ Pd- Shown are the results for 
Pl = 75, pD = 25 {Pl/Pd = 3) in blue, pL = 100, Pd = 25 {Pl/Pd = 4) in 
black, and Pl = 125, Pd = 25 [Pl/Pd = 5) in red. In this case the net chiral 
excess at steady-state is ecss — 0.38, eess — 0.47, and ecss — 0.53, respectively 
(corresponding to 38%, 47%, and 53% enantiomeric excess). The choice of a is 
predicated by the value of Pl{d), since for a <^ Pl(d) the suppression of activated 
monomers would be too effective. We investigated values for a = 100 and a = 300 
but did not find any significant variation in the outcomes. 

These test runs show that for differing left and right reaction rates a substantial 
amount of chiral excess may be reached at steady-state. Although for these exam- 
ples the differences were fairly large, we note that if the chiral-selective changes 
in rates appear in exponential factors (as in the example above for temperature 
dependence), small changes in the parameters may generate fairly large changes 
in the resulting reaction rates. 



• u/. = 5,ao = 10, 


= ^5o = 100 


• ai_ = S.ao = 15, 


= = 100 


• at ^ 5,(^0 ^20, ISl 








= aD = 


10, /?i=75, Ad = 25 


• at 




10, ^i. = 100. = 25 


• ffL 




10, pL - 125, l3o - 25 




time (d s ^) 



time (ds ^) 



Fig. 1 Time evolution of enantiomeric excess for illustrative examples of reaction rates with 
cr = 200. Left: (Si^ = /Sjj = 100, with varying ratios of a^Ja^. Shown are = 5.0, ajj = 10.0 
(a^/ao = |) in blue, = 5.0, q_d = 15.0 = |) in black, and = 5.0, a^, = 20.0 

{oil/<^d = j) in r^d. Right: = cto = 10.0, with varying ratios of Pd/I^l- Shown are 
= 75, /3_D = 25 (liL/Po = 3) in blue, Pl = 100, = 25 [Pl/Pd = 4) in black, and 
Pl = 125, Pd = 25 [Pl/Pd = 5) in red 



In general, we note the reactions favor an L-excess for ratios az, /aD < 1 and/or 
Pl/Pd > 1 as demonstrated by the results shown in Figure 1. At first glance this 
may seem counterintuitive, particularly the relationship between oll/ old and the 
net enantiomeric excess of homochiral polymers. Intuitively one would expect that 
increasing ud (the polymerization rate of activated D-monomers) while holding 
OLL (the polymerization rate of activated L-monomers) fixed would result in a 
net excess favoring D, and not L as observed. However, on closer inspection one 
must note that old describes the rate of attachment of activated _D-monomers to 
growing _D— or L- chains, leading to higher yields of _D-enriched mixed polymers 
and a correspondingly lower yield of pure D-polymers. Likewise, increasing the 
deactivation rate of L-monomers (i.e. Pl^ and thereby increasing the ratio Pl/Pd 
leads to fewer L-activated monomers for polymerization which decreases the yield 
of pure homochiral L- polymers, but also decreases the loss of L-monomers to 
polymers of mixed chirality, resulting in a net increase in the enantiomeric excess 
favoring L. Clearly, with the same values of L and D reaction rates, ee = for all 
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t. We must thus examine how the resuhs behave for statistical samplings of a^, 
an, Pl, and Pd- 

It is also important to compare the production of homochiral and heterochiral 
polymers, since the naive expectation is that the system will be swamped by 
heterochiral polymers. In Table [l] we used the same reaction rates as in Figure [l] 
to show that this is not the case. In fact, the ratio of homochiral to heterochiral 
polymers [defined as 7?, = {C + T)) / Ad] is fairly invariant with respect to changes in 
the reaction rate parameters (grouped in the Table as blocks of three rows each), 
suggesting the possible existence of attractors in the phase space of the reaction 
kinetics determining the equilibrium concentration of polymers. 





do 


Pl 


Pd 


ee 


M 


C + V 


n 


5 


10 


100 


100 


0.242 


49.154 


8.328 


0.170 


5 


15 


100 


100 


0.384 


47.365 


7.255 


0.153 


5 


20 


100 


100 


0.479 


46.345 


6.630 


0.143 


10 


10 


75 


25 


0.376 


33.958 


3.391 


0.100 


10 


10 


100 


25 


0.466 


36.743 


3.931 


0.107 


10 


10 


125 


25 


0.530 


39.105 


4.405 


0.113 



Table 1 Comparison of homochiral to heterochiral polymer concentrations in steady state. 
The top three rows correspond to the parameters of the left plot in Figure [l] while the bot- 
tom three rows corresponds to the parameters of the right plot. We introduced the ratio of 
homochiral to heterochiral polymers, TZ = {C + 'D)/JV{. For each of the two groups (top three 
and bottom rows) note how the ratio is quite insensitive to variations in the reaction rates. 



We have performed a detailed statistical study of the reaction equations with 
fluctuating values for the L- and D- reaction rates q_l, qd, Pl and Pd as described 
in eqs. [5] and [6] scaled by d. This involves solving the coupled system of ordinary 



differential equations above (eqns. 24 30) N times, each with values for the four 



reaction rates given hy a.L = a -\- 5l, an = a -\- 5d, Pl = P + and P = 
Pd + (,D, where the bars denote the mean values (cf. eq. [s]) and 5i^(^d) and ^l(D) 
are Gaussian-distributed random numbers (cf . eq. [6| , within a fixed rms width set 
by Qo and Pq. 

Our results indicate that the net chirality is overall more sensitive to the ampli- 
tude of the chiral-selective variations about the mean values of the reaction rates 
than to the values of the reaction rates themselves. In Figures 2 and 3 we show the 
results for three sets of data with variations in Pl(d) between 10% and 30% (i.e. 
0.1/3 < ^0 < 0.30/3) and in Ql(u) fixed at 20% (i.e. ao = 0.2a). Each data point 
represents an ensemble average over a random sample of iV = 5000 individual 
experiments (or prebiotic realizations). Shown are results for a = 10, /3 = 100 in 
blue, a = 15, /3 = 90 in red, and a = 20, /3 = 120 in black. Figure [2] demonstrates 
that the ensemble-averaged enantiomeric excess is racemic, with an equal number 
of systems favoring excesses in D as in L, as should be expected. Closer inspec- 
tion of the details of the distribution demonstrate that some chemistries within 
the statistical sampling have moderate to large enantiomeric excesses. In Figure 
3a, the mean value of the net chirality is shown (i.e. the mean magnitude of the 
enantiomeric excess). Although the distribution of enantiomeric excess is centered 
on zero as shown in Figure [2j Figure 3a demonstrates that the average net chi- 
ral asymmetry obtains moderate values, with the (|ee|) varying from ~ 10% for 
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/3o = 0.1/3 to ~ 16% for f3o = 0.3^, with the (|ee|) increasing nearly linearly with 
the size of the fluctuations. We stress that individual members of the ensemble 
may have much larger net chiralities than these mean values. The variances of 
the data set are shown in Figure 3b, and range between ~ 0.015 and ~ 0.045 for 
po = 0.1/3 and /3o = 0-313, respectively. As with the mean value of the net chirality, 
the variance increases with the amplitude of departure from the mean, leading to a 
wider spread in the values of the enantiomeric excess with increasing /?o . For large 
departures from the mean, say /3o = 0.3/3, the variance is ~ 0.04, corresponding 
to a standard deviation of 0.2, indicating that 68% of systems have an |ee[ < 0.2, 
and 27% have an enantiomeric excesses in the range 0.2 < |ee| < 0.4 (assuming a 
normal distribution for the data, see Figure 4). A small fraction of 0.2% of systems 
possess a sizable |ee| in excess of 0.6. 

In Figure [4j we further explore the spread in the distribution of the enan- 
tiomeric excess in experimental systems for an ensemble with a = 20, j3 = 120 
and ao = 0.25q! and /3o = 0.25/?. Outliers in the distribution have very high enan- 
tiomeric excesses up to as much as 80 — 90%, although most systems fall within a 
range of enantiomeric excesses with |ee| < 0.25, where the mean of the distribution 
lies. In Figurejsjwe repeat the analysis with a = 20, P = 100 and ao = O.lOa and 
/3o = 0.10/3. Note that even with rms of 10%, it is still possible to obtain substan- 
tial enantiomeric excesses above 20% at a few percent. Some outliers obtain over 
40% ee. 

In Figure [6] we show a surface plot of the mean values of the net chiral asym- 
metry for I3l{D) varying between 10% and 30% (i.e. 0.1/3 < /3o < 0.3/3) and 
fluctuations in ciLi^D) varying over the same range. The highest values for the net 
enantiomeric excess are observed for large fluctuations in a^j-jj) and /?£,(£)), where 
mean values for the net enantiomeric excess reach as high as |ee| ~ 0.2 for the sys- 
tems under study. These results demonstrate a general trend of increased average 
(|ee|) with increased chiral selectivity in the rates, as well as an increase in the 
spread of the distribution, indicating that larger departures form the mean lead 
to a larger fraction of outliers with significant enantiomeric excess. 
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Fig. 2 Ensemble averaged enantiomeric excess, (ee), for sample sizes = 5000 and fluctu- 
ations in ctj^^^ij fixed at 20% (ao = 0.2(5), with fluctuations in Pl(D) varying between 10% 
and 30% (0.1/3 < /3o < 0.3/3) and fixed cr = 200. Shown are « = 10, ^ = 100 in blue, a = 15, 
/3 = 90 in red, and a = 20, /3 = 120 in black. 
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Fig. 3 Ensemble averaged mean of the net enantiomeric excess, |ee|, and variance for the 
ensemble data shown in Figure j2j with A'^ = 5000, with fluctuations in ct^^jj) fixed at 20% 
(ao = 0.2(5), and fluctuations in P^i^jj) varying between 10% and 30% (O.l^S < /3o < 0.3;S) 
and fixed a = 200. Shown are a = 10, /3 = 100 in blue, a = 15, ;3 = 90 in red, and a = 20, 
/3 = 120 in black. Left: Mean values of the net enantiomeric excess, {|ee|). Right: Variance 
in the enantiomeric excess ((ee — (ee))^). 
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Fig. 4 Distribution of ensemble ee values for A'' = 5000 with a = 20, /3 = 120, «o = 0.25q: 
and /3o = 0.25^, and cr = 200 . 



4 Summary and Outlook 

Due to the symmetry between the chemistries of L- and D- monomers, including 
the biologically relevant examples of amino acids and nucleotides, it is expected 
that nearly all chemical interactions between the two enantiomeric forms and their 
environment are symmetric with respect to both chiralities. However, specific lo- 
cal environmental influences, such as mineral surfaces, chiral-selective temperature 
sensitivity, photolysis, or the presence of asymmetric catalysts, can bias reactions 
to favor one chirality over the other to a larger or lesser extent. All of these exam- 
ples are inherently globally symmetric, such that on average a statistical sampling 
of chemistries will be racemic (i.e. some will favor L and others D). Asymmetry 
therefore arises depending on the particular context. For example, surfaces can 
lead to chiral-selectivity due to differential activity of the two enantiomers when 



adhered to minerals or clays (Hazen and ShoU 2003 Joshi, Pitsch, and Ferris 
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Fig. 5 Distribution of ensemble ee values for TV = 1 X lO'' with a = 20, /3 = 100, ao = 0.10a, 
/3o = 0.10/3, and cr = 200 . 
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Fig. 6 Ensemble averaged net enantiomeric excess, (|ee|), for N = 100 with fluctuations in 
'^L{D) 3.nd Pl{D) between 10% and 30% (ao/a, Po/P between 0.1 and 0.3) with a = 20, 
B = 200 and a = 200. 



2007). This can lead specific local chemistries to have a chiral bias in reaction 



rates that differ significantly from the standard reactivities in free solution. We 
have explored general scenarios of a similar nature here, within the context of the 
as yet unknown prebiotic chemistry which may have given rise to the first life- 
forms. As we have argued, different environmental contexts will lead to varying 
deviations from the mean reactivity of L- and D- enantiomers, leading to local- 
ized symmetry breaking. Larger fluctuations will result in some chemistries having 
substantial chiral asymmetries. Based on the results presented here, we see that 
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fluctuations as smaU as 10% in reactivities of L- and D- monomers may lead to 
significant chiral symmetry breaking {i.e. |ee| > 10%). 

Additionally, fluctuations in chiral-selective reaction rates may provide an al- 
ternative mechanism to the initiation of chiral asymmetry prior to ampliflcation 
via autocatalysis with enantiomeric cross-inhibition in a Frank-type mechanism. 
Traditional approaches assume that fluctuations in the relative numbers of L- and 
D- monomers are enough to break the initial symmetry between the two enan- 
tiomers, which may then b e amplified. However, fluctu ations in chemical systems 
typically scale as 1/ \/\N) (Landau and Lifshitz 19801, where A'' is the number of 
molecules in the system, and will therefore be exceedingly small for most chemical 
systems. Deviations from perfect racemates will therefore be correspondingly small 
and require significant enzymatic enhancement by efficient catalysts to be ampli- 
fied {e.g. see Gleiser and Walker (2008|). Here we present an alternative and more 



robust mechanism for symmetry breaking, whereby the initial symmetry break- 
ing occurs at a macroscopic rather than microscopic level, being a manifestation 
of fluctuations in reactivities, rather than in molecular numbers. Such scenarios 
obviate the need for special asymmetric initial conditions, allowing asymmetry to 
emerge as the reactions unfold, due to asymmetries present in the local chemistry. 

Finally, we note that our results suggest that the origin of homochirality on 
the prebiotic Earth was a statistical process, potentially hosting many abiogenic 
events, few of which were highly asymmetric. Quite likely, only those chemistries 



which were highly asymmetric were capable of giving rise to life ( jBonner 1995 



Fitz et al. 2007), and only one of these trials ultimately led to the emergence 



of the LUCA. This has different observational consequences for experimental sys- 
tems as well as searches for biosignatures on other worlds than the exogenous or 
endogenous sources of chiral excess described above. Parity-violation would result 
in a universal bias toward a particular chiral asymmetry, while UV-irradiation 
might lead to chiral biases on the scale of planetary systems. Here, chiral selec- 
tivity is a local statistical phenomenon, and we should expect a large sampling of 
extraterrestrial stereochemistries to be racemic on average. 
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